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Abstract 

(-H We prove a constant term theorem which is useful for finding weight 

polynomials for Ballot/Motzkin paths in a strip with a fixed number 

S of arbitrary 'decorated' weights as well as an arbitrary 'background' 
weight. Our CT theorem, like Viennot's lattice path theorem from 
which it is derived primarily by a change of variable lemma, is ex- 
i— I pressed in terms of orthogonal polynomials which in our applications 

of interest often turn out to be non-classical. Hence we also present 
an efficient method for finding explicit closed form polynomial expres- 
sions for these non-classical orthogonal polynomials. Our method for 
finding the closed form polynomial expressions relies on simple com- 
binatorial manipulations of Viennot's diagrammatic representation for 
orthogonal polynomials. In the course of the paper we also provide 
Q\ a new proof of Viennot's original orthogonal polynomial lattice path 

theorem. The new proof is of interest because it uses diagonalization 
of the transfer matrix, but gets around difficulties that have arisen in 
• *h past attempts to use this approach. In particular we show how to sum 

over a set of implicitly defined zeros of a given orthogonal polynomial, 
cither by using properties of residues or by using partial fractions. We 
conclude by applying the method to two lattice path problems impor- 
tant in the study of polymer physics as models of steric stabilization 
and sensitized flocculation. 

Keywords: Lattice Path, Dyck Path, Ballot Path, Motzkin Path, 
Paving, Three-term Recurrence, Jacobi Matrix, Tri-diagonal Matrix, Trans- 
fer Matrix, Orthogonal Polynomial, Chebyshev Polynomial, Rational Gen- 
erating Function, Constant Term, Rogers formula, Asymmetric simple ex- 
clusion process. 
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1 Introduction and Definitions 



As is well known to mathematical physicists the form of the solution to a 
problem is often more important than its existence. Such is the case in this 
paper. Determining the generating function of Motzkin path weight polyno- 
mials in a strip was solved by a theorem due to Viennot [29] (see Theorem 
[T] below - hereafter referred to as Viennot 's Theorem). In the applications 
discussed below what is required are the weight polynomials themselves. 
Whilst these can be written as a Cauchy integral of the generating function 
this form of the solution is of little direct use for our applications. To this 
end we have derived a related form of the generating function given by The- 
orem [2] (hereafter referred to as the Constant Term, or CT, theorem). The 
CT theorem is certainly well suited to extracting the weight polynomials 
for Chebyshev type problems (section [7| and as a starting point for their 
asymptotic analysis |25j . 

The CT theorem is proved, as we shall show below, by starting with Vi- 
ennot 's Theorem and using a 'change of variable' lemma (see Lemma 5.1). 
We will also provide a new proof of Viennot 's Theorem that is based on 
diagonalizing the associated Motzkin path transfer matrix. The latter proof 
is included as it rather naturally leads to the CT Theorem. It also has some 
additional interest as it has several combinatorial connections [7[. For ex- 
ample, a combinatorial interpretation of what has previously appeared only 
as a change of variable to eliminate a square root in Chebyshev polynomials 
turns out to be the generating function of binomial paths. Another connec- 
tion is a combinatorial interpretation of the Bethe Ansatz [1] as determining 
the signed set of path involutions, as for example, in the involution of Gessel 
and Viennot [20 in the many path extension [8 . 

Two classes of applications for which the CT theorem is certainly suited 
are the Asymmetric Simple Exclusion Process (ASEP) and directed mod- 
els of polymers interacting with surfaces. For the ASEP the problem of 
computing the stationary state, and hence in finding phase diagrams for 
associated simple traffic models, can be cast as a lattice path problem 
[21 El IH El HU1 HH HS1 US]. For the the ASEP model the path problem 
required is actually a half plane model with two weights associated with the 
lower wall (the upper wall is sent to infinity to obtain the half plane.) |10j . 
In chemistry the lattice paths are used to model polymers in solution |15] - 
for instance in the analysis of steric stabilization and sensitized flocculation 

CUE]. 

In the application section of this paper we find explicit expressions for 
the partition function - or weight polynomial - of the DiMazio and Rubin 
polymer model [T7] . This model was first posed in 1971 and has Boltzmann 
weights associated with upper and lower walls of a strip containing a path. 
The wall weights model the interaction of the polymer with the surface. 
The solution given in this paper is an improvement on that published in 
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|7J. Previous results on the DiMazio and Rubin model have only dealt with 
special cases of weight values, for example, the case where a relationship 
exists between the Boltzmann weights [S] . We also present a natural gener- 
alisation of the DiMazio and Rubin weighting which may have application 
to models of polymers interacting with colloids [23], where the interaction 
strength depends on the proximity to the colloid. 

In order to use the CT Theorem both the ASEP and polymer models 
require explicit expressions for 'perturbed' Chebyshev orthogonal polynomi- 
als. Their computation is addressed by our third theorem, Theorem [3] 

1.1 Definitions and Viennot's Theorem 

Consider length t lattice paths, p = VQViV2---Vt, in a height L strip with 
vertices Vi G S = Z>o X {0, 1, L}, such that the edges e» := Vi-iVi satisfy 
V{ — Vi-i G {(1, 1), (1, 0), (1, —1)}. An edge is called up if Vi — = (1, 1) 
, across if Vi — Vi-\ = (1,0) and down if Vi — V{-\ = (1,-1). A vertex 
Vi = (x, y) has height y. Weight the edges according to the height of their 
left vertex using 



w(ei) 



1 if ei is an up edge 

6fc if ei is an across edge with Vi-i = {i — 1, k) (1) 
Afc if ei is a down edge with Vi_i = (i — 1, k). 



The weight of the path, w(p), is defined to be the product of the weights of 
the edges, ie. for path p = voVi...vt, 

t 

w(p) = Y[w(ei). (2) 

8=1 

Such weighted paths are then enumerated according to their length with 
weight polynomial defined by 

Z t (y',y;L):=Y,w(p) (3) 



where the sum is over all paths of length t, confined within the strip of height 
L, y' is the height of the initial vertex of the path and y is the height of the 
final vertex of the path. An example of such a path is shown in Figure [T] 
These paths are weighted and confined elaboration's of Dyck paths, Ballot 
paths and Motzkin paths - for enumeration results on these classic paths 
see, for example, [22J. Once non-constant weights are added, many classical 
techniques do not (obviously) apply. 

This work focuses upon solving the enumeration problem for the types of 
weighting in which a small number of weights take on distinguished values 
called 'decorations', and the rest of the edges have constant 'background 
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Figure 1: An example of a lattice path of length 15, in a strip of height 
L = 2, with weight ^o^i^l^i-^i anc ^ starting at y' = and ending at y = 1. 



weights' of just one of two kinds, accordingly as the step is an across step 
or a down step. 

We introduce notation to describe the positions of the decorated edges: 
let 2?b C {0, L] and V\ C {1, L} be sets of integers called respectively 
decorated across-step heights and decorated down-step heights. 
Then paths are weighted as in Equation (Jl|, with 



b ifi(£V b 

b + bi \iieV b 

A ifi^A 

\ + if i G V x 



(4a) 
(4b) 



where b and A will be called background weights and b\ and \ will be 
called decorations. 

The generating function for the weight polynomials is given in terms 
of orthogonal polynomials by a theorem due to Viennot j29l 130] - see also 

nam] 



Theorem 1 ( [2H] ) • The generating function of the weight polynomial ^ 
given by 



IS 



M L (y',y;x) := £ Z t (y', y; L) x l 



„Y-Y 



,R Y >(x) h y ,, y R^\x) 



Rl+i(x) 



(5) 



where h y > iy = H y<l < y > A; if y' > y and h y ,, y = 1 if y' <y,Y' = min{y',y} 
and Y = max{y', y}. The polynomials Rk{x) are the reciprocal polynomials 
and Ru\x) the shifted reciprocal polynomials defined by 



R k (x) := x k P k {l/x) 



and 



R { k j \x) := R k (x) 



i+j 



(6) 



Xi — > x 



where the orthogonal polynomials P k {x) satisfy the standard three term re- 
currence jig Pffi . 



Pk+i{x) = {x-b k )P k {x)-\ k P k -i{x), P (x) = 1, Pi(x) = x-b Q . (7) 
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This theorem may be proved in several ways; as the ratio of determi- 
nants - see |27j section 4.7.2, by continued fractions |19| I2T] or by heaps of 
monomers and dimers [3T] . In Section [4] we provide a new proof that uses 
diagonalization of the transfer matrix of the Motzkin paths. 

2 A Constant Term Theorem 

Our main result is stated as a constant term of a particular Laurent expan- 
sion. Since the constant term method studied in this paper depends strongly 
on the choice of Laurent expansion, we briefly recall a few simple facts about 
Laurent expansions. Since we only consider rational functions we restrict 
our discussion to them. A Laurent expansion of a rational function about a 
point z = Zi is of the form ^2 n>no a n (z — Zi) n . The coefficients a n depends 
on the chosen point z% and the annulus of convergence. Furthermore, the 
nature of no generally depends on three factors: i) no > (ie. the series is a 
Taylor series) if the annulus contains no singular points, ii) no < is finite 
if the inner circle of the annulus contains only a non-essential singularity at 
z% and iii) no = — oo if the inner circle contains at least one other singularity 
at z Zi or an essential singularity. 

In this paper we only need case ii), with z% = and the series convergent 
in the annulus closest to the origin. Thus the constant term is defined as 
follows. 

Definition 1. Let f(z) be a complex valued function with Laurent expansion 
of the form 

oo 

/(*) = a nZ n (8) 

n=riQ 

with no S Z. Then the constant term in z of f{z) is 

CT z [f(z)] = oo. (9) 

This is, of course, just the residue of f{z)/z at z = 0. Note, the form 
of the Laurent expansion given in Q uniquely specifies that it corresponds 
to that Laurent expansion of f{z) that converges in the innermost annulus 
that is centred at the origin. 

Our main result gives the weight polynomial for Motzkin paths in a 
strip as the constant term of a rational function constructed from Laurent 
polynomials. The Laurent polynomials we use, (p), are defined in terms 

of the conventional (shifted) orthogonal polynomials, pjf\x), by the simple 
substitution 

Lf{ P ):=P^{x{p)) (10) 

with 

x(p) = p + b + \p-\ (11) 
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The orthogonal polynomials Pk(x) = P[ (%) satisfy the standard three term 
recurrence Q which, for the shifted polynomials Pjf\x), becomes 

pW( x ) = ( x _ 6fc+i _ 1 )p fc W( a; ) - A fc+j _iP^ 2 (x), fc > 2 (12) 



Pi 



P (J) (^) 



X 

1 



and Afe ^ V/c. We now state our principal theorem. 

Theorem 2 (Constant Term). Let Zt{y' ,y; L) be the weight polynomial for 
the set of Motzkin paths with initial height y' , final height y, confined in a 
strip of height L, and weighted as specified in Equations Q. Then 



P 



Z t (y',y;L) = CT p [p + b + 
with Y' = minjy', y},Y = max{y', y} 
h 



A\* L Y '{p) hy>, y L^y\p) 



"y ,y 



Tly<l<y' 

1 



Ll+i(p) 



ify'>y 

otherwise. 



(13) 



(14) 



(i) 

and the Laurent polynomials LY given by (10) 



The form of this constant term expression should be carefully compared 
with that arising from Viennot's theorem when used in conjunction with the 
standard Cauchy constant term form for the weight polynomial (for y' < y), 



Z t (y',y,L) = CT X 



1 



r-i+l 



X 



y _ yl R y ,(x)Rty \x) 
Rl+i(x) 



(15) 



In particular, (13) is not obtained by simply substituting 1/x = p + b + Xp 



into (15) as was done to define the Laurent polynomials L 



It is a useful exercise to compare the difference in effort in computing a 
general expression for Z t (y',y;L) in the simplest possible case y = y' = 0, 
bfr = 1, Afc = A starting from (|13l) compared with starting from (15). 



Assuming that a simple explicit expression is desired for the weight poly- 
nomial then the utility of the CT Theorem depends on and arises from three 
factors. The first problem is how to calculate the orthogonal polynomials. 
For various choices of the weights and A^ , the classical orthogonal polyno- 
mials are obtained and hence this problem has already been solved. However 
for the applications mentioned earlier, which require 'decorated' weights, the 
polynomial^] do not fall into any of the classical classes. Thus computing 
the polynomials becomes a problem in itself and is addressed by Theorem |3j 

lr They may be thought of as 'perturbed' Chebyshev polynomials 
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The second problem is more subtle and is concerned with how the poly- 
nomials are represented. Whilst Pk(%) is by construction a polynomial this 
is not necessarily how it is first represented. For example, consider a Cheby- 
shev type polynomial, which satisfies the constant coefficient recurrence re- 
lation 

Sk+i(x) = xSk(x) - Sk-i(x), Sx(x) = x, So(x) = l. (16) 

This recurrence is easily solved by substituting the usual trial solution Sk = 
v k with v a constant, leading immediately to a solution in the form 

af . (x + Vx^4) k+1 - (x - Vx^4) k+l , 1V . 
2 k+1 Vx 2 - 4 

Whilst this is a polynomial in x, in this form it is not explicitly a polynomial 
as it is written in terms of the branches of an algebraic function. The repre- 
sentation of the polynomials is important as it strongly influences the third 
problem, that of computing the constant term (or residue). If the polyno- 
mials are explicitly polynomials (rather than, say, represented by algebraic 
functions) then the obvious way of computing the Laurent expansion, and 
hence residue, is via a geometric expansion of the denominator. Whilst in 
principal this can always be done for a rational function the simpler the 
denominator polynomials the simpler the weight polynomial expression - in 
particular we would like as few summands as possible, preferably a number 
that does not depend on L, the height of the strip. 

The fact that this can be achieved for the applications studied here shows 
the advantage of the CT Theorem in this context over, say, the Rogers 
Formula [26] (see Proposition 3A of [19]), in which the weight polynomial is 
always expressed as an L-fold sum no matter how small the set of decorated 
weights. 



3 A Paving Theorem 

Our second theorem is used to find explicit expressions for the orthogonal 
polynomials which are useful to our applications. These are polynomials 
arising from problems where the number of decorated weights is fixed (ie. 
independent of L). Theorem [3] of which we make extensive use, expresses 
the orthogonal polynomial of the decorated weight problem in terms of the 
orthogonal polynomials of the problem with no decorated weights (ie. Cheby- 
shev type polynomials). 

Theorem 3 (Paving). 1. For each c S {1, 2, k — 1}, we have an 'edge 
cutting' identity 

P^(x) = P^\x)Pl + :\x) - X c+j P^P^ix) (18a) 
2 The square roots can of course be Taylor expanded to show explicitly it is a polynomial. 
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and a 'vertex-cutting' identity, 



pjf\x) = (x- b c+j ) P^(x)Pjf_ + c - + 1 1 \x) 

- \ j+c+1 P^\x)P^:+ 2 2 \x) - X c+j P?\{x)P^\x), (18b) 

where P^\x) satisfies ( |l~2"| ). 

2. Fix j and k. Let \T>b\ and \T>\\ be the number of decorated 'across' and 
'down' steps respectively whose indices are strictly between j — 1 and 
j + k. Let d = \V X \ + \ V b \. Then 

3 max Imax 

Pi 3 \x) = ^a J l[S k ^x), (19) 



3=1 i=l 



where 



l<jmax<2 lVx] Z lVb] , l<W<d + l; (20) 

and kji is a positive integer valued function; decorations are all con- 
tained in the coefficient dj's, and the S^^x)^ are the background 
weight dependent Chebyshev orthogonal polynomials satisfying 

Sk+i{x) = {x-b) S k {x) -XSk-i(x), (21) 

with Si(x) = x — b, So(x) = 1 and A / 0. 

We do not give an explicit expression for the kji as it is strongly depen- 
dent on the sets T>i, and T>\. They are however simple to compute in any 
particular case, for example see (58) in the application section [71 The sig- 
nificance of (19) is that it shows that the decorated polynomials P^ can be 
explicitly expressed in terms of the undecorated (ie. Chebyshev) polynomials 

The first part of the paving Theorem [3] follows immediately from an 
'edge-cutting' and a 'vertex-cutting' technique respectively, applied to Vi- 
ennot's paving representation of orthogonal polynomials, which we describe 
in Section [7] This geometric way of visualising an entire recurrence in one 
picture is powerful; from it we see Part [T] of the Theorem as a gestalt, so in 
practice do not need to remember the algebraic expressions but may work 
with paving diagrams directly. Part [2] follows immediately by induction on 
Part[T] 



4 Proof of Viennot's theorem by Transfer Matrix 
Diagonalization 

In this section we state a new proof of Viennot's theorem. This proof starts 
with the transfer matrix for the Motzkin path (see section 4.7 of [27] for 



8 



an explanation of the transfer matrix method) and proceeds by diagonal- 
izing the matrix. As is well known this requires summing an expression 
over all the eigenvalues of the matrix. The eigenvalue sum is a sum over 
the zeros of a particular orthogonal polynomial. This sum can be done for 
the most general orthogonal polynomial even though the zeros are not ex- 
plicitly known. We do this in two ways, the first uses two classical results 



(Lemma 4.2 and Lemma 4.3) and a paving polynomial identity. The essen- 
tial idea is to replace the sum by a sum over residues and this residue sum 
can then be replaced by a single residue at infinity. The second proof uses 
partial fractions. 

The transfer matrix, Tl, for paths in a strip of height L, is a square 
matrix of order L + 1 such that the (y',y) th entry of the t th power of the 
matrix gives the weight polynomial for paths of length t, i.e. 

Z t (y',y;L) = (Ti) y ,, y . (22) 

For Motzkin paths with weights Q the transfer matrix is the Jacobi matrix, 



(bo 


1 











o\ 


Ai 


h 


1 

















h 


1 




















bh-2 


1 

















b L -i 


1 


\o 











Al 


bj 



(23) 



The standard path length generating function for such paths, with spec- 
ified initial height y' and final height y, is given in terms of powers of the 
transfer matrix as 



M L (y', y; x) := £ Z t (y', y; L) x l = £(T|; 



y ,y 



(24) 



which is convergent for \x\ smaller than the reciprocal of the absolute value 
of the largest eigenvalue of Tl- 

The details of the proof are as follows. We evaluate Z t {y' , y; L) 
by diagonalization 

T l L = VD l L U (25) 

with Dl = diag(xo, x±, ■■■,xl) a diagonal matrix of eigenvalues of Tl, and 
V and U respectively matrices of right eigenvectors as columns and left 
eigenvectors as rows normalized such that 



UV = I, 



(26) 



where I is the unit matrix. One may check that this diagonalization is 
achieved by setting the i th column of V equal to the transpose of 



v(xj-i) = (Po(xi- 1 ),P 1 (Xi-l},...,P L (xi-.l)) 



(27) 
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and the i th row of U equal to 



u(xi-i) 



X\...\l 



' Poixi-x) P^Xi-l) P 2 {Xi-i) 



P L (Xj-l) 

1 Ai A1A2 ' Ai...Al 

where the set of eigenvalues {xi}f =0 are determined by 

P L+1 {xi) = 0. 



(28) 



(29) 



with the orthogonal polynomial Pl+i given by the three term recurrence 



(12). Orthogonality of left with right eigenvectors of the Jacobi matrix (23) 



follows by using the Christoffel-Darboux theorem for orthogonal polynomi- 



als. Equation (26) then follows as (29) gives L + 1 distinct zeros and hence 



L + l distinct eigenvalues and hence L + l linearly independent eigenvectors. 

For simplicity in the following we only consider the case y' < y in which 
case the h y ^ y factor is one - it is readily inserted for the case y' > y. Thus, 
multiplying out Equation (25) and extracting the (y',y) th entry, we have 



( T l)y',y - (Aj,+i-..Al) 



X \Pyl (Xj)Py{Xi) 

^ P' L+l {xi)P L ( Xi ) 



(30) 



Note that P' L+1 (xi) 7^ and Pl(xi) 7^ by the Interlacing Theorem for 
orthogonal polynomials, so that all the terms in the sum are finite. Since 
Viennot's theorem does not have a product of polynomials in the denom- 



inator we need to simplify (30), which is achieved by using the following 
lemma. 



Lemma 4.1. Let Xi be a zero of Pl + \{x). Then 

X y+1 ...X L Py(Xi) = PUx^P^iXi) 



(31) 



This lemma follows directly from the edge-cutting identity (18a) by 
choosing k = L + 1, j = c and c = L — h together with the assumption 
that Pl+i(xi) = 0, to obtain a family of identities parametrized by h; which 
are then iterated with h = 0, 1, . . . ,L — y — 1. 



Applying Lemma 4.1 to (30) gives the following basic expression for the 



weight polynomial resulting from the transfer matrix 



( T l)y',y 



E 



(32) 



Note, this use of the transfer matrix leads first to an expression for the 
weight polynomial, to get to Viennot's Theorem we still need to generate on 
the path length and also simplify the sum over zeros. The former problem is 
trivial, the latter not. We do the sum over zeros in two ways, first by using 
a contour integral representation and secondly by using partial fractions. 
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Using a contour integral representation 



To eliminate the derivative in the denominator of ( 32 ) we use the following 
Lemma. 

Lemma 4.2. Let P{z) and Q{z) be polynomials in a complex variable z and 
P(zi) = 0, then 



Res 



P'(z) 
Q(*)p^,{*,*} 



(33) 



where mi is the multiplicity of the root Z{, and P'{z) is the derivative with 
respect to z. 

The following Lemma allows us to replace a residue sum with a residue 
at infinity. 

Lemma 4.3. Let P{z) and Q{z) be polynomials in a complex variable z, 
then 



1 



2m L Q(z) 



dz 



^2 Res 

Res 



Q(z) 

POM 
z 2 Q(l/z) 



{z,Zi} 
{z,0} 



(34a) 



(34b) 



where 7 is a simple closed anticlockwise-oriented contour enclosing all the 
zeros of Q(z) and A is the set of zeros ofQ(z). 



Note, (|34j) simply states that the sum of all residues of a rational function, 
including that at infinity, is zero. These lemmas are proved in most books 
on complex variables - see j5] for example. 

Now use Lemma |4.2| to get rid of the derivative in the denominator 
(rrtj = 1 as all zeros of orthogonal polynomials are simple), to produce 



(Tt) 



{xi\P L+1 (xi)=0} 



Res 



(</+!), 



X 



Pl+i(x) 



{x, x{\ 



(35) 



Applying Lemma |4.3| to sum over the zeros, gives the weight polynomial as 
a single residue (or constant term) 



Dy',y 



Res 



P y ,(l/ X )Pt + y\l/x) 

X*+ 2 P L+1 (1/X) 



{x,0} 



(36) 



As noted above, a factor, h y ^ y needs to be inserted in the numerator for the 
case y 1 > y. Comparing this form with ( 15 ) and changing to the reciprocal 
polynomials ^ gives Viennot's theorem. Thus we see the change to the re- 
ciprocal 1/x in this context corresponds to switching from a sum of residues 
to a single residue at infinity. 
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Using partial fractions 

We can also derive a residue or constant term expression for the weight 



polynomial, (36), without invoking either of the calculus Lemmas 4.2 



or 



Lemma 4.3 by using a partial fraction expansion of a rational function. In 
particular, if we have the rational function 



G(x) 



Q{x) 

T( x y 



with Q and T polynomials of degree a and b respectively and a < b. Assum- 
ing T(x) is monic and has simple zeros T{xi) = we have T(x) = Yli(x-Xi) 
and thus we have, using standard methods, the partial fraction expansion 



i 



where T'{x) is the derivative of T(x) and thus 

i 

Geometric expanding each term gives us the coefficient of x n in G(l/x) as 

[x n ]G(l/x)=^Xr 1 §^y 



If we now compare (32) with (38) we see that Q — ► Py'Pij—y ■> T 
and t = n — 1, thus we get 

,(»+!)/ 



(38) 



M P y/ (l/x)P [ L y l 1 '(l/x) 



Res 



i p,,(i/x)pt + ; ) (i/x) 

.'+ 2 Pl+i(1/x) 



{x,0} 



(39) 



Note, the orthogonal polynomials satisfy the conditions required for the 



existence of (37). Thus we see that the sum over the zeros (ie. eigenvalues) 



in (32) is actually a term in the geometric expansion, or Taylor series, of a 



partial fraction expansion and hence 'summed' by reverting the expansion 
back to the rational function it arose from. 

Note, although the partial fraction route is elementary, to get from the 
natural representation of the t th power of a matrix in terms of its diagonal- 
ization (ie. equations (25) and (|30|)) to Viennot's theorem we still needed 
Lemma 14.11 
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5 Proof of the Constant Term Theorem 



For the proof of the CT Theorem we use the following residue 'change of 
variable' lemma. 

Lemma 5.1 (Residue change of variable) . Let f{x) andr(x) be a functions 
which have Laurent series about the origin and the Laurent series of r{x) 
has the property that r(x) = x k g{x) with k > and g(x) has a Taylor series 
d( x ) = En>o a nX n such that clq / 0. Then 



Res[f(x),x] = ~^^ es 



(40) 



where Res[f(x),x] denotes the residue of f{x) at x = 0. 



The lemma appears to have first been proved by Jacobi |23j. A proof 
may also be found in Goulden and Jackson |21| . Theorem 1.2.2. Note, the 
condition on r(x) is equivalent to the requirement that the Laurent series of 
the inverse function of r{x) exist. 

To prove the CT Theorem we start with Viennot's theorem and use the 
simple fact that the coefficient of x l in the Taylor series for f{x) is given by 



CT S 



1 



x 



Res 



1 



rt+l 



/(*) 



(41) 



which gives (15), that is, 



Z t (y',y;L) = CT X 



1 y _ y , Ryi(x)h yl , y R>l_ 



(y+i), 



x 



t+l' 



Rl+i(x) 



Now consider the change of variable defined by 

x{p) = 2 P , . , 
p z + bp + A 

which has the Taylor expansion about the origin 

,-,2 



(42) 



(43) 



x(p) 



P _ 

A A 2 



and thus x(p) satisfies the conditions of Lemma |5,l| with k = 1. Note, (43) 
is the reciprocal of the change given in ([ll]). Thus from (42) and using 
Lemma 5.1 with the change of variable (|43[) and derivative 



dp 



x{p) 



x P - 



1 



(p + 6 + Ap-i) 



1A2- 



(44) 



we get the CT theorem result (13), in terms of the Laurent polynomials 
defined by dlob. 
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We make two remarks. The first is the primary reason for the change of 
variable in this instance is that it "gets rid of the square roots" such as those 



that appear in the representation (17) and hence changes a representation of 



a Chebyshev type polynomial in terms of algebraic functions to an explicit 
Laurent polynomial form. 

For the second remark we note that for this proof of the CT Theorem 
we could equally well have started just before the end of the diagonalization 



proof of Viennot theorem, ie. equation (36), and proceeded with the residue 



change of variable (43). In other words, the transfer matrix diagonalization 
naturally ends with an expression for the weight polynomial - this has to 
be generated on before getting to Viennot 's theorem, whilst the first step in 



this CT proof, ie. equation (42), is to undo this generating step. 



6 Proof of the Paving Theorem 

The paving interpretation of the orthogonal polynomial three term recur- 
rence relation was introduced by Viennot |29j. We will use this interpreta- 
tion as the primary means of proving Theorem [3j First we define several 
terms associated with pavings. A path graph is any graph isomorphic to 
a graph with vertex set {vi}*_ and edge set {^TO+iI^Tq 1 . A monomer is 
a distinguished vertex in a graph. A dimer is a distinguished edge (pair 
of adjacent vertices). A non-covered vertex is a vertex which occurs in 
neither a monomer nor a dimer. A paver is any of the three possibilities: 
monomer, dimer or non-covered vertex. A paving is a collection of pavers 
on a path graph such that no two pavers share a vertex. We say that a 
paving is order k if it occurs on a path graph with k vertices. An example 
of a paving of order ten is shown in Figure [2] 



I I I 



Figure 2: The upper part of this diagram shows a paving. The lower part 
indicates the reason for calling it a 'paving' , as it is in bijection with a more 
standard 'paving' or 'tiling' diagram, where long and short tiles are used, 
the short tiles being of two possible colours. 

Weighted pavings are pavings with weights associated with each paver. 
We will need pavers with shifted indices in order to calculate the shifted 
paving polynomials that occur in Theorem [3] Thus, the weight of a paver 
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a with shift j is defined as follows. 

uij(a) - 



x the paver a is a non-covered vertex 

—bi+j the paver a is the monomer (45) 
— the paver a is the dimer ej 



The weight of a paving is defined to be the product of the weights of the 
pavers that comprise it, i.e. 

for p a paving. It is useful to distinguish two kinds of paving. Pavings 
containing only non-covered vertices and dimers are called Ballot pavings; 
those also containing monomers are called Motzkin pavings. A paving 
set is the collection of all pavings (of either Ballot or Motzkin type) on a 
path graph of given size. We write 

Vk &l = {p\p is a Ballot Paving of order k} (47) 
-pMotz = £p| p ig a Mo t z kin Paving of order k} (48) 

When it is clear by context whether we refer to sets of Ballot or Motzkin 
pavings, the explanatory superscript is omitted. 

A paving polynomial is a sum over weighted pavings defined by 



with Vk being either "PP al or "Pj^ otz and weights as in Equation ( 45 ) . If 



k 



k = then we define Pq (x) := 1 Vj. The diagrammatic representation for 
a paving polynomial on a paving set of Ballot type is 

pj?\x)<- c>^±^4^o cy-^¥^ty-^=b , Kn , 

"0 »1 "2 Vk-2 (5U) 

where the question mark denotes that the edge can be either a dimer or not. 
The diagrammatic representation for a paving polynomial on a paving set 
of Motzkin type is 

P [ k 3 \x)^ 0-^H3M^2) <!y^¥^y-=^b , K1 v 

"0 Oi «2 Ofe_2 Vk-1 ("''-J 

Note, we overload the notation for pf?\x) since we will use it for the set 
of pavings and the corresponding paving polynomial obtained by summing 
over all weighted pavings in the set. 



Viennot [29] has shown that Equation (49) satisfies (12) and hence 



pjf\x) is an orthogonal polynomial. Ballot Pavings correspond to the case 
bi = Vi. An example of Ballot pavings and the resulting polynomial is 
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-> X A - (Ai + A 2 + X 3 ) x 2 + AiA 3 
and an example of a Motzkin paving set with associated polynomial is 

p 3 (6) (x) <- ^fe4 8 = 




-> x 3 - (b 6 + b 7 + b 8 )x 2 

-(A 7 + A 8 - &6&7 - b 6 b s - b 7 b 8 )x 
-(b 6 b 7 b 8 - 6 6 A 8 - 6 8 A 7 ) 

The two identities of Theorem 2 correspond to 

• cutting at an arbitrary edge; and 

• cutting at an arbitrary vertex. 

In the first procedure we consider the edge e c . This edge is either paved 



or not paved with a dimer. Equation (52) illustrates the division into these 



two cases, and the corresponding polynomial identity. Note that the right 
hand side of the expression obtained is a sum of products of smaller order 
polynomials such that the weight '— A c+ j' associated with edge e c occurs 
explicitly as a coefficient; and is not hidden inside any of the smaller order 
polynomials. 



>u) 



Pl J '(x) <- ® © ?^2> T^-CD f^O) ® 

= ® ® — ?— ^— <V) © — ^±^±^2; © 



^(^^(x) - \ c+j P?\{x)Pt^\x) (52) 
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The second procedure is to cut at an arbitrary vertex 'i> c '. In this procedure 
the cases to consider are: v c is non-covered, v c is a monomer, v c is the 
leftmost vertex of a dimer, and v c is the rightmost vertex of a monomer. 



These four cases are shown in Equation (53), and the resulting identity 



gives Pu (x) as a sum of four terms, each of which contains a product of 
two smaller order polynomials. The vertex weight '— 6 C +/ occurs explicitly 
as a coefficient, and is not hidden in any of the smaller order polynomials. 



-bc+j-2 _^ —l>c+j-i_^ —b c +j _^ —bc+j+i 

— © — ¥^kD — ¥^<D — f^y- 



-bj+k-i 

-CD 



= ®- 

-bj 
+ ®- 



+ 



+ 



~bj 

©- 



j 



-Oc+j-2 _\ —Oc+j-1 

— © ^=te 

-bc+j-2 _ x -bo+j-i 

— © -f^D 

-bc+j-2 \ —Oc+j-1 

~~ ac+i— 



X -bc+j+1 

o ©— 



Oc+j+1 




-6j+fc-l 



°c+j+2-0j + k-l 



-3 



b c+j )P<?Hx)PtJ-x\x) X c+j+1 pV\x)Ptl?\x) 

-K +j P^\{x)P^\x) 



(53) 



These two 'cutting' procedures, as shown in Diagrams (52) and (53), prove 



Part [T] of Theorem [3j Finally Part [2] of the Theorem follows immediately 
from Part [T] by induction on the number of decorations of each of the two 
possible kinds ('across step' and 'down step'). 

We note that the proof of Theorem [3] shows that the upper bounds given 
in Equations (20) are tight only when the decorations are well-separated 



from each other as well as from the ends of the diagram. Otherwise, pulling 
out a given decoration as a coefficient may pull out a neighboring decoration 
in the same procedure, meaning that fewer terms are needed. Thus we 
can deal more efficiently with decorated weightings in which collections of 
decorations are bunched together, as illustrated in the examples in Section [Tj 



7 Applications 

We now consider two applications. The first is the DiMazio and Rubin 
problem discussed in the introduction. Solving this model corresponds to 
determining the Ballot path weight polynomials with just one upper and one 
lower decorated weight. The second application is an extension of that prob- 
lem in which two upper and two lower edges now carry decorated weights. 
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7.1 Two decorated weights 



For the DiMazio and Rubin problem we need to compute the Ballot path 



weight polynomials with weights given by ( 54 ) . A partial solution was given 



in [12] for the case when the upper weight is a particular function of the 
lower one ie. k + u = kuj. The first general solution was published in [7] in 
2006. The solution we now give is an improvement on [7j (which was based 




#— « 



Figure 3: An example of a two weight Dyck path (above) in a strip of height 
three and the corresponding paving problem (below). 



on a precursor to this method) as it contains fewer summations. 



Theorem 4. Let Z2 r (n,uj; L) be the weight polynomial for Dyck paths of 
length 2r confined to a strip of height L, and with weights (see Qj, b = A = 
1, bi = \/i and 



uJ 
K ■ 




ifi = L 
ifi = \ 
otherwise. 



(54) 



Then 



where 



Z2r(K, UJ] L) = CT 



(P + P 



-l\2r, 



Ap L - Bp 



-L 



ACp L - BDp- 



(55) 



A 
B 
C 
D 



p — UJ 
I -Cop 2 

2 



P~ 
1 



Kp 



(56a) 
(56b) 
(56c) 
(56d) 



k := k — 1 and Co := uj — 1. 
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An example of the paths referred to in Theorem [4] is given in Figure [3j 
Theorem [4] is derived using Theorem [2] and Theorem |3j We shall show the 
derivation of the denominator of ( |55[ ) but omit the details for the numerator 
which may be similarly derived. 



We work with polynomials in the x variable, and then use Equation ( 10 ) 



to get the expression in terms of p. First write Pl+i{x) for our given weight- 



ing as in Equation (57): 



(57) 



We cut at the first and last edges, as in Equation (58). (The edges to 



cut at are chosen since they mark the boundary between decorated and 
undecorated sections of the path graph.) 



Pl+i(x) 



-l -l -l 



-l -l -l 



L-4 L-3 L-2 L-l L 




o 



o 



L-4 L-3 L-2 L-l L 



-1 -1 -1 



-1 -1 




Now we have Pl+i(x) as a sum of four terms, as in Equation ( |59[ ): 
Pl+i{x) = x 2 S L -i(x) - kxS L - 2 {x) - ojxS L - 2 (x) + kojS L -3(x). 



(58) 



(59) 



Each of the four terms is a product of three contributions - the X S. K s 
and u's come from the short sections in each row of the diagram, and S^s 
represent the long sections, as in Equation (60): 

Sk(x) = ° °~~^9 9 Pr^a^S^i (60) 



Now Equation (60) represents polynomials satisfying the recurrence and 
initial conditions given in Equation (21 ), so we may substitute Equation (17) 
for each occurrence on an 'SV in Equation (59). We now have a sum of 
ratios of surds, which may be simplified by the change of variables specified 
in Equation (10) to yield 

R L+1 (p) =P L+1 (p + p- 1 ) (61) 
= {ACp L -BDp- L )p 2 {p-p- 1 ) (62) 
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which is, up to a factor, the denominator of (55 ). The numerator is similarly 
derived. 



Next we indicate how to expand Equation ( 55 ) to give an expansion for 



the weight polynomial in terms of binomials. Our particular expression is 
obtained via straightforward application of geometric series and binomial 
expansions, with the resulting formula containing a 5-fold sum. It is cer- 
tainly possible to do worse than this and obtain more sums by making less 
judicious choices of representations while carrying out the expansion, but it 
seems unlikely that elementary methods can yield a smaller than 5-fold sum 
for this problem. 



First, manipulate the fractional part of Equation (55) into a form in 



which geometric expansion of the denominator is natural, as in Line (63) 
below. Do the geometric expansion and multiply out to give two terms: 
extract the m = case from the first term and shift the index of summation 



in the second to give Line ( 65 ) . 



Ap L - Bp~ L 
ACp L - BDp- L 



D 



bdP 



AC J2L 

bdP 



(63) 



1 

D 



A 

B~T) i 



„2L 



E 

m=0 



2mL 



(64) 



D 



+ E 



m=l 



P 



2mL 



E 

m=l 



mr*im— 1 



A m C 



J2mL 



B m D r 



(65) 



We next find the CT, separately, of each of the three terms multiplied by 
series for (p + p~ ) 2r (l — p 2 )- The first gives 
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CT 



(p + p-i)*-(l-/) 



(66) 



CT 



' 2r fir"- 
\u=0 



(1 " P 2 ) 



\ ~~ ~m „2m 



\m=0 



K p 



(67) 



J2 * m CT 



m=0 



2r 

E 

u=0 



2r 



(p 



2r-2u+2m _ ^2r-2u+2m+2j 



(68) 



5> 



m=0 



2r 
r + m 



2r 

r + m + 1 



(69) 



m=0 

(2n\ I 2n 



(70) 



where C n -k := (u) — (&— i)' The secon d term is expanded similarly, with 
the positive powers of A and C and the negative powers of B and D each 
contributing a single sum, which, when concatenated with the original sum 
over m, creates a 5-fold sum altogether. The third term generates a similar 
5-fold sum; and this difference of a pair of 5-fold sums is combined into one 
in the final expression in Theorem [5] below. 

Theorem 5. Let Z2 r (n,uj;L) be as in Theorem^ The 



len 



m>0 



7/1 

^ _ E < 

m>l pi,P2>0 si,S2=0 



]^Sl+S2£S2+P2 ( 2 ) Sl+Pl 



-> E 

>P2> 

m\ ( m\ (m — 1 + pi \ (m + p2 



(71) 



si/ \s 2 

C r -r—k—l 



Pi 

m — S2 
m+p 2 



P2 



C r 



;i — k 



where k = p\ +P2 — si — S2 + (L + 2)m — 1, and C n; fc zs the extended Catalan 
number, 

CTH-. 

T7ie binomial coefficient (") is assumed to vanish if n < or m < or 
n < m. 
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Note, when trying to rearrange ( 71 ) care should be taken when using 



any binomial identities because of the vanishing condition on the binomial 
coefficients - the support of any new expression must be the same as the 
support before (alternatively the upper limits of all of the summations must 
be precisely stated). 

7.2 Four Decorated weights 

This second problem is a natural generalization of the previous problem. We 
now have a pair of decorated weights in the pair of rows adjacent to each 
wall, as in Figure [4] In the earlier DiMazio- Rubin problem, paths have been 
interpreted as polymers zig-zagging between comparatively large colloidal 
particles (large enough to be approximated by fiat walls above and below) 
with an interaction occurring only upon contact between the surface and 
the polymer; this weighting scheme could be used to model such polymer 
systems, but now with a longer range interaction strength that varies sharply 
with separation from the colloid. 



/ 



V 



u>2 



\ 



• 

«1 



Figure 4: An example of a four weight Dyck path (above) in a strip of height 
five and the corresponding paving problem (below). 



Theorem 6. Let Zi t {k\, K2> w%\ L) be the weight polynomial for Dyck 
paths of length 2r confined to a strip of height L, and with weights (see Q 
fe = A = 1, 6i = 0Vi and 



( 

U>l 


- 1 


ifi = L 


Ll>2 


- 1 


ifi = L- 


< K2 


- 1 


ifi = 2 


Kl 


- 1 


ifi = \ 




V 




otherwise. 



(72) 
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Then 



where 



Z 2r (Kl,K2,0Jl,U2;L) 



CT 



(p + p 



A 
A 
B 
B 

C 

c 



-l\2r 



f ABp L — A Bp 



-L 



{ 



CBp L -CBp~ 



(P' l -P) 



1 
1 

P 



K2P 
K-2P 2 

(&l +L02)P~ 



- U2P 

P^ 1 - (0>l + UJ2)p - U)2P 3 

p-(k x + k 2 )p~ 1 - k 2 p~' 
p' 1 - («i + k 2 )p - k 2 p 3 



for ki 



1 and Coi 



1. 



(73) 



(74a) 
(74b) 
(74c) 
(74d) 
(74e) 
(74f) 



An example of the paths referred to in the above theorem is given in 
Figure |4} 

The constant term expression in Theorem[6]may be expanded in a similar 
fashion to that of Theorem|4]to yield a 9-fold sum. The fractional component 
of Equation ( 73 ) may be written 



ABp L - ABp 
CBp L -CBp- 



-L 



A 



m=l 



C m+1 B'' 



2mL 



m=l 



2mL 



(75) 



by a method precisely analogous to that applied in the 2-weights case. When 
multiplied by {p+ p~ 1 ) 21 {p^ 1 — p) and the constant term extracted, the initial 



term gives the double summation in the first term of Equation ( 76 ) in The- 



orem [7] below. The other two terms of Equation (75) each give 9-fold sums, 
as a consequence of the double sums yielded by each of the powers, positive 
and negative, of C, B, C and B . These two 9-fold sums are combined in 
the second term of Equation ( 76 ) in Theorem [7] 
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Theorem 7. Let Z 2r (K\, k%,ui,u<2,; L) be as in Theorem^ Then 
[ 

i>0 j=0 



"0 + 2/ V "o + 1/ \u 



m s\ m S2 

' Si \ / m \ / S2 \ I v\ 



EEEEEEEEE 

m>l si=0 ii=0 s 2 =0 i 2 =0 Vl >Q ji=0 v 2 >0 j 2 =0 

v 2 + m-l\ fv 2 _ +S2+n+i2>< 



m - 1 / VJ2 

K 2 (Ki + K2) U! 2 \U)\ -\- U)2) x 

'";0(^: m ^ 1 )(^, 2 - 1 )- fe+l >(:) + U 1 



(76) 



for uq = r + 2i — j and u\ = r + mL + v\ + V2 + s\ + S2 +ji +j 2 — 2i\ —1i2- 

Theorems [5] and [7] are to be compared with the Rogers formula below, 
which gives the weight polynomial as an order L-fold sum. 

Theorem 8 (Rogers |26j). Let Z2n be the weight polynomial for the set of 
Dyck paths of length 2n with general down step weighting (see Qj.' b = A = 
1, bi = Vi and Aj = Kj — 1 in either a strip of height L or in the half plane 
(take L = oo). Then the weight polynomial is given by 

min{n— 1,L— 1} 

Z 2n (K 1 ,K2,...;L) = ^2 si, (77) 

1=0 

where si is the weight polynomial for that subset of paths in Z2 n which reach 
but do not exceed height L + 1. The s; 's are given by 

so = «?, (78) 

and 

jo-l Ji-l Ji-i-l l-l , . . 1X 

* = £ E ••• E nn -T yr h 4~ h ---^ l+ \ (79) 

h=l 32=1-1 31=1 k=0 V Jfe J 
for I > 1; with 

jo ■= n, (80) 
ji+i ■= 0. (81) 
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We see that there is a trade-off between having comparatively few sums 
(compared with the width of the strip) but a complicated summand, as in 
Theorems [5] and [7j versus having a simpler summand but the order of L 
sums irrespective of the number of decorations. 

It would be an interesting piece of further research to see whether the 
solutions to the problems presented in Section [7J containing as they do mul- 
tiple alternating sums, are in some appropriate sense best possible or not. 
Another potentially useful area of further research would be an investiga- 
tion of good techniques for extracting asymptotic information directly from 
the CTp expression. It would also be interesting to have a pure algebraic 
formulation of this constant term method - one that does not rely on any 
residue theorems. 

Acknowledgements 

Financial support from the Australian Research Council is gratefully ac- 
knowledged. One of the authors, JO thanks the Graduate School of The 
University of Melbourne for an Australian Postgraduate Award and The 
Centre of Excellence for the Mathematics and Statistics of Complex Sys- 
tems (MASCOS) for additional financial support. We would also very much 
like to thank Ira Gessel for some useful insights and in particular to drawing 
our attention to the Jacobi and Goulden and Jackson reference [22] used in 
the proof in Section [5j 

References 

[1] H. A. Bethe. Z. Phys., 71:205, 1931. 

[2] R A Blythe and M R Evans. Nonequilibrium steady states of matrix- 
product form: a solver's guide. J. Phys. A: Math. Theor, 40:R333-R441, 
2007. 

[3] R. A. Blythe, M. R. Evans, F. Colaiori, and F. H. L. Essler. Exact 
solution of a partially asymmetric exclusion model using a deformed 
oscillator algebra. arXiv:cond-mat/9910242, 2000. 

[4] R. A. Blythe, W. Janke, D. A. Johnston, and R. Kenna. Dyck paths, 
motzkin paths and traffic jams. arXiv:cond-mat/0405314, 2000. 

[5] Ralph P. Boas. Invitation to complex analysis. Random House, New 
York, 1987. 

[6] R. Brak, S. Corteel, J. Essam, R. Parviainen, and A. Rechnitzer. A 
combinatorial derivation of the PASEP stationary state. Electronic 
Journal of Combinatorics, (R108), 2006. 



25 



[7] R. Brak, J. Essam, J. Osborn, A. Owczarek, and A. Rechnitzer. Lattice 
paths and the constant term. Journal of Physics: Conference Series, 
42:47-58, 2006. 

[8] R. Brak, J. Essam, and A. L. Owczarek. From the bethe ansatz to the 
gessel- viennot theorem. Annals of Comb., 3:251-263, 1998. 

[9] R. Brak, J. Essam, and A. L. Owczarek. Exact solution of N directed 
non-intersecting walks interacting with one or two boundaries. J. Phys. 
A., 32(16):2921-2929, 1999. 

[10] R Brak and J W Essam. Asymmetric exclusion model and weighted 
lattice paths. J. Phys. A: Math. Gen., 37:4183-4217, 2004. 

[11] R Brak, A L Owczarek, and A Rechnitzer. Exact solutions of lattice 
polymer models. J. Chem. Phys., 2008. Submitted. 

[12] R Brak, A L Owczarek, A Rechnitzer, and S G Whittington. A directed 
walk model of a long chain polymer in a slit with attractive walls. J. 
Phys. A: Math. Gen., 38:4309-4325, 2005. 

[13] T. S. Chihara. An introduction to orthogonal polynomials, volume 13 
of Mathematics and its Applications. Gordon and Breach, 1978. 

[14] S. Corteel, R. Brak, A. Rechnitzer, and J. Essam. A combinatorial 
derivation of the PASEP algebra. In FP 'SAC 2005. Formal Power Series 
and Algebraic Combinatorics, 2005. 

[15] P.-G. de Gennes. Scaling Concepts in Polymer Physics. Cornell Uni- 
versity Press, Ithaca, 1979. 

[16] B Derrida, M Evans, V Hakin, and V Pasquier. Exact solution of a Id 
asymmetric exclusion model using a matrix formulation. J. Phys. A: 
Math. Gen., 26:1493 - 1517, 1993. 

[17] E A DiMarzio and R J Rubin. Adsorption of a chain polymer between 
two plates. J. Chem. Phys., 55:4318-4336, 1971. 

[18] M. R. Evans, N. Rajewsky, and E. R. Speer. Exact solution of a cellular 
automaton for traffic. J. Stat. Phys., 95:45-96, 1999. 

[19] P. Flajolet. Combinatorial aspects of continued fractions. Discrete 
Math., 32:125-161, 1980. 

[20] I. M. Gessel and X. Viennot. Determinants, paths, and plane partitions. 
Preprint., 1989. 

[21] I. P. Goulden and D. M. Jackson. Path generating functions and con- 
tinued fractions. J. Comb. Th. Ser. A, 41:1 - 10, 1986. 



26 



[22] P. Goulden and D. M. Jackson. Combinatorial Enumeration. John 
Wiley & Sons Inc, 1983. 

[23] C.G.J. Jacobi. De resolutione aequationum per series infinitas. Journal 
fur die reine und angewandte Mathematik, 6:257-286, 1830. 

[24] A L Owczarek, R Brak, and A Rechnitzer. Self-avoiding walks in slits 
and slabs with interactive walls. J. Chem. Phys., 2008. Submitted. 

[25] A. L. Owczarek, T. Prellberg, and A. Rechnitzer. Finite-size scaling 
functions for directed polymers confined between attracting walls. J. 
Phys. A: Math. Theor., 41:1-16, 2008. 

[26] L. J. Rogers. On the representation of certain asymptotic series as 
continued fractions. Proc. Lond. Math. Soc.(2), 4:72-89, 1907. 

[27] R Stanley. Enumerative Combinatorics: Vol 1, volume 1 of Cambridge 
Studies in Advanced Mathematics 49. Cambridge University Press, 
1997. 

[28] G. Szego. Orthogonal Polynomials. Amer. Math. Soc, 1975. 

[29] G. Viennot. A combinatorial theory for general orthogonal polynomials 
with extensions and applications. Lecture notes in Math, 1171:139-157, 
1985. 

[30] X. Viennot. Une theory combinatoire des polynomes orthog- 
onaux generaux. Unpublished, but may be downloaded at 
http://web.mac.com/xgviennot/Xavier Viennot /livres. html, 1984. 

[31] X. G. Viennot. Heaps of pieces, i: Basic definitions and combinatorial 
lemmas. Lecture notes in Math, 1234:321, 1986. 



27 



